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Three-particle azimuthal correlation measurements with a high transverse momentum trigger par- 
ticle are reported for pp, d-l-Au, and Au-I-Au collisions at ^s„^ = 200 GcV by the STAR experiment. 
Dijet structures are observed in pp, d+Au and peripheral Au-|- Au collisions. An additional structure 
is observed in central Au-|-Au data, signaling conical emission of correlated charged hadrons. The 
conical emission angle is found to be ^=1.37 ± 0.02(stat)^g Q^(syst), independent of p±. 
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Collisions at BNL's Relativistic Heavy Ion Collider 
(RHIC) create a hot and dense medium that cannot 
be described by hadronic degrees of freedom l]. Evi- 
dence of this is provided, in part, by jet-quenching: on 
the away side of a high transverse momentum {p±_) trig- 
ger particle (in azimuth relative to the trigger particle, 
A(j) = (j) — (f>t ~ tt) , the correlated yield is strongly sup- 
pressed at > 2 GeV/c 2], while at lower p_\_ the 
yield is enhanced and the correlated hadrons appear to 
be partially equilibrated with the bulk medium and are 
broadly distributed in azimuth 3] . A number of physics 
mechanisms may account for the data: broadened jets 
due to large angle gluon radiation deflected jets due 
to collective radial flow of the bulk [Bl or pathlength 
dependent energy loss and conical emission due to 
Cerenkov gluon radiation 0] or Mach-cone shock waves 
generated by large energy deposition in the hydrody- 
namic medium 

Identifying the underlying mechanism is important as 
it may probe the medium properties such as its speed 
of sound and equation of state @, 0|. To discriminate 
between the various mechanisms, we have performed an 
analysis of three-particle azimuthal correlations between 
a high trigger particle and two lower p±^ associated 
particles in A0i = 0i — 0t (*=lj2) We integrate over 
the pseudo-rapidity {r]) direction because the near- and 
away-side jets are not correlated in 77 Many mech- 
anisms predict that pairs of associated hadrons will be 
shifted away from A^= tt, but will remain close to each 
other (A(?!)i « A^a) In contrast, the Mach-cone 

or Cerenkov radiation scenarios would result in particle 
emission on a cone around the away-side jet axis. When 
projected onto the azimuthal direction, the strongest sig- 
nal of conical emission would be Jacobian peaks where 
pairs of correlated hadrons appear with equal probabil- 
ity to be close together or to be far apart and symmetric 
about TT (i.e., A(?!)i - tt w tt - A(/)2) % El. The lat- 
ter feature is specific to conical emission. In this let- 
ter, we present evidence for correlated hadron pairs that 
are symmetrically located about tt relative to the trig- 
ger particle. The analysis is carried out with trigger 
and associated particles of 3 < pj_ < 4 GeV/c and 
1 < Pj_ < 2 GeV/c, respectively, in pp, d-|-Au, and 
Au+Au colhsions at y/s^ — 200 GeV. 

Details of the STAR (Solenoidal Tracker at RHIC) ex- 
periment are described elsewhere 12]. This analysis uses 
2x10^ pp, 6.5x10^ d-fAu, and 1.2x10'' minimum bias 
(MB) and 1.9x10^ central trigger Au-|-Au events taken 
in 2001-2004. The central trigger data set corresponds to 



approximately 12% of the total geometric cross-section, 
and will be henceforth referred to as "12% central" colli- 
sions. Charged particles are reconstructed with the Time 
Projection Chamber (TPC) [13|, which sits in a uniform 
0.5 Tesla magnetic field. The Au-|-Au data are divided 
into nine centrality bins according to the charged parti- 
cle multiplicity in the pseudorapidity region |7y|<0.5 as 
in [l4i- Similarly the d-t-Au data are divided into three 
centrality bins of 0-10%, 10-20%, and 20-100%. The 
trigger and associated particles are restricted to |''?|<1. 
Our results are corrected for the centrality-, pi_-, and 
0-dependent reconstruction efhciency for associated par- 
ticles and the 0-dependent efhciency for trigger particles, 
and are normalized per corrected trigger particle. 

Various approaches r nay be taken to measure three- 
particle correlations [13, El, [3] . This analysis treats the 
event as composed of two components: one is correlated 
with the trigger, I2, and the other is background uncor- 
related with the trigger except the indirect correlation 
via anisotropic flow. The correlated particle distribution 
(two-particle correlation) is thus given by 



y2(A0) = r2(A0) - aBinci^2(A(/.) , 



(1) 



where l2(A(/)) = dN/dAcf) is the raw associated particle 
density per trigger. The other, background term is con- 
structed by mixing triggers with different inclusive events 
(i.e. MB events within a given centrality bin), with the 
effect of anisotropic flow, 

F2(A(/)) = 1 + 2v!^*\2 cos(2A0) + 2vf\i cos(4A0) , (2) 

constructed pair-wise using flow measurements (vn^ and 
Vn, n = 2,4, are trigger and associated particle n**^ har- 
monic coefficients, respectively) [l^ 12 1; Anc = Mnc/27r 



is the inclusive event associated multiplicity density; 
a = Nhg/Ninc scales iVinc to the underlying background 
associated multiplicity A'^bg in trigger events, as discussed 
below. 

In our two- component approach, the full three-particle 
distribution, ¥3 , consists of the correlated triplets of in- 
terest, I3, sets of three particles that are uncorrelated 
with each other except via flow, and cases where two 
of the particles are correlated (including jets and other 
correlations such as resonance decays) and the third is 
uncorrelated with the first two except via flow. The cor- 
related pair distribution (three-particle correlation) is ob- 
tained by [13, [nl 



F3(A01,A,/)2) = y3(A'^l,A02)-aBinc[^'2(A0i)F2(A02)+r2(A,^2)^^2(A0i) 

F3(A0i,A02^ 



-6a^r2'"'(A<^i,A02) 



F2(A</)i - A(^2) 



(3) 
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where Y3 = cP N / dA(j)idA(j)2 is the raw associated par- 
ticle pair density per trigger, and the second and third 
terms on the r.h.s. are backgrounds. The second term, re- 
ferred to as trig-corr-bkgd, arises from combining a corre- 
lated trigger-associated pair with a background particle, 
and is constructed from the product of the two-particle 
correlation and its flow modulated background. 

The third term, referred to as trig-bkgd-bkgd, arises 
from combining a trigger with two background particles, 
and contains all correlations between the two background 
particles as well as the flow correlation between them and 
the trigger. The former is the inclusive event pair density 
1^2'"'^ = d^N[nc/dA(j)idA(j)2 relative to a random trigger 
0t, which is constructed by mixing the trigger from one 
event with two particles from another, inclusive event. 
The latter is referred to as trigger flow, where 

FsiAcbuAcj)^) = f^2(A0i) + J^2(A02)-2 

+ 2v^^\^^\^^^ cos2(A0i - 2A02) 



+ 2vX'v. 



(*)„(2)„,(i) 



cos2(2A0i - A02) 



+ 2i;^^^i;^^^4*^cos2(A0i -f A02) (4) 

is constructed triplet-wise by mixing the trigger with 
particles from two different inclusive events. The fac- 
tor ba? scales the number of pairs in inclusive events, 
(Mnc(Mnc — 1)), to that in the underlying background, 
{Nhg{N\)g — 1)). Non-Poisson multiplicity fluctuations, 
which can be different in inclusive events and in the back- 
ground underlying trigger events, result in deviations of 

b from one. We approximate b by ^jv- ^^{N^~-iy)f{N- — y 
where N is the associated multiplicity in trigger events. 

The analysis procedure is performed and the scaling 
factors a and b are determined for each centrality bin sep- 
arately; the final three-particle results are combined over 
centrality bins to increase the statistics. The value of a is 
determined assuming that the three-particle correlation 
signal has Zero Yield at Minimum (3-ZYAM); the total 
size of the minimum signal regions is chosen to be 10% 
of (27r)^. It is so chosen so that it is small enough to ap- 
proximate the real minimum, but large enough to avoid 
large statistical fluctuations. This size is varied between 
5-15% of (27r)2, keeping a fixed, to assess the systematic 
uncertainty on b. The upper end of the systematic uncer- 
tainty on a is taken to be the a value from two-particle 
ZYAl (Zero Yield At 1 radian) where l2(A0) vanishes at 
|A^±1| <7r/18 0. The lower end is determined, while 
keeping b at its default value, from the lowest data point 
(out of total 24x24), which should be lower than the 
true 3-ZYAM because of statistical fiuctuations. With 
a at each systematic end, the value of b is readjusted, 
shifting the three-particle correlation result by an ap- 
proximately constant pedestal, to preserve 3-ZYAM. For 
the top 5% centrality fraction with the 12% central data, 
a = 0.994(irool) and b = 1. 00021 (t°:°°g[|^) ^ 

Figure [T] shows two-particle correlations in Au-|-Au 



central collisions: the raw Y2{A(j)) and the a-scaled back- 
ground aBinc^2(A(/)) in (a), and the background sub- 
tracted Y2{A(j)) in (b). Fitting l2(A^) to various func- 
tional forms similar to those in Ref. [3] yields away- 
side peaks centered 1.18-1.34 radians from tt. Fig- 
ure[ljc,d,e) show, respectively, the raw three-particle cor- 
relation A(j)t , A(j>2) , fca^K)'"'^, and the trig-corr-bkgd 
term plus trigger flow |19j . 

Table U summarizes the major sources of systematic 
uncertainties. (I) Uncertainty in the normalization fac- 
tor a is assessed as above. (II) The V2 used is the av- 
erage of modifled reaction plane W2{mrp} and 4-particle 
cumulant ■i;2{4} 0|. The two-particle cumulant V2{2}, 
which contains flow fluctuations and potentially non-flow 
effects, and the i;2{4} or V2{2d} (obtained from a two- 
dimensional analysis in A77 and A0) bracket the sys- 
tematic uncertainties. We used a parameterization of 
V4 = 1.15 u| [l3| and the V2 uncertainties are propa- 
gated. (Ill) The trig-corr-bkgd term in Eq.® is con- 
structed from the two-particle correlation and its back- 
ground, both averaged over the reaction plane (RP) an- 
gle. The effect of the change of the correlation structure 
with the angle between the trigger and the RP [i^] is 
estimated and included in our final results. The size of 
the effect is assigned as a single-sided systematic uncer- 
tainty. The systematic uncertainty from (I) primarily 
impacts the overall magnitude of the correlation, with 
little infiuence on the shape, whereas those from (II) and 
(III) have a smaller impact on the magnitude, but affect 
the shape of the correlation. 

Table HI also lists the total systematic uncertainty from 
other, minor sources: uncertainty in the normalization 
factor 6 estimated as above; ±20% uncertainty on the 
uncertainties due to the finite cen- 



unmeasured 



(t) 



17| 



trality bins on trig-corr-bkgd and trig-bkgd-bkgd terms 
estimated by breaking each centrality into finer bins; and 
10% uncertainty due to the efficiency correction. 



TABLE L Systematic uncertainties on three-particle correla- 
tion strength on the away side: central region (|A(jf>i,2 — 7r| < 
0.35) and ofT-diagonal region (|A(^i — tt ± 1.37| < 0.35 and 
|A<^2 - TT =F 1-371 < 0.35). 



source 


d-l-Au 


50-30% Au-hAu 


12% central Au-I-Au 




cent. 


cent. 


off-diag. 


cent. 


ofT-diag. 


(I) a 


+160/ 

-18 /O 


+29 (V 
-60 /o 


+30 (V 
-63/0 


+42 (V 
-61 /o 


+210/ 


(II) V2 




+ 17% 


+36 0/ 


-130/ 


+32 (V 
-16'° 


(III) RP 




+ 11% 


+8% 


+32% 


+5% 


others 


+ 110/ 
-10 


+ 120/ 
_ii/o 


+ 180/ 
-12/0 


+160/ 


±12% 



Figure [5] shows the background subtracted three- 
particle azimuthal correlations. Is, in MB pp, d±Au, and 
three combined centrality bins of MB Au±Au and the 
12% central collisions. Four distinct peaks are observed 
for each data set, corresponding to both correlated par- 
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FIG. 1: (a) Raw two-particle correlation signal Y2 (red), background aBincF2 (solid histogram), and background systematic 
uncertainty from a (dashed histograms), (b) Background-subtracted two-particle correlation Y2 (red), and systematic uncer- 
tainties due to a (dashed histograms) and flow (blue histograms), (c) Raw three-particle correlation Y3. (d) ba^Y2'^'^. (e) Sum 
of trig-corr-bkgd and trigger flow. Data are from 12% central Au-|-Au collisions. Statistical errors in (a,b) are smaller than the 
point size. 




Ad) =(|) -d) 

^1 ^1 ^Trig 



FIG. 2: Background subtracted three-particle correlations, 
%, for (a) pp, (b) d+Au, (c) 80-50% Au-^Au, (d) 50-30% 
Au+Au, (e) 30-f0% Au-FAu, and (f) 12% central Au-I-Au. 
Statistical errors per bin are approximately ±0.012 in (a) and 
±0.006 in (b), both at (tt.tt), and are ±0.022, ±0.049, ±0.099 
and ±0.077 from (c) to (f), similar for all bins. 



tides on the near side (around A(j>i — A02 = 0), both 
on the away side (around tt), and one on each side. The 
near-side peaks are slightly elongated along the diago- 
nal, probably due to momentum balance in combination 
with the fact that the trigger direction differs from its 
parent's. 

The away-side central peak is elongated along the di- 
agonal, progressively from pp to d+Au to Au+Au colli- 



sions. This indicates that the away-side pairs stay rel- 
atively close while their angles vary over a wide range 
event-by-event. Figure [3l^a) shows the effect quantita- 
tively by projecting the d+Au three-particle correlation 
on the away side (1 < A(/)i.2 < 27r — 1) along the di- 
agonal in S = (A(/)i + A(j>2)/2 — tt and off-diagonal in 
A = (A01 - A(/)2)/2, within the ranges of < A < 0.35 
and < 0.35, respectively [l9[. For comparison the 
off-diagonal projection on the near side (|A0i.2| < 1) is 
also shown. 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 



-2 -1.5 -1 -0.5 0.5 1 1.5 2 
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FIG. 3: Projections of away-side three-particle correlations 
along the diagonal E within < A < 0.35 (squares) and along 
the off-diagonal A within |Ej < 0.35 (points) in (a) d+Au and 
(b) 12% central Au+Au collisions. The shaded areas indicate 
systematic uncertainties on the off-diagonal projections. The 
histogram in (a) is the near-side off-diagonal projection. The 
histogram in (b) is the away-side off-diagonal projection of 
our result with a = 6 = 1. 

For central Au+Au collisions, additional peaks are ob- 
served in Fig. [2] on the away side along the off-diagonal, 
indicating large opening angles between the away-side 
correlated pairs, symmetric about tt, A(/)i — vr « tt — A(j)2 
corresponding to each off-diagonal peak. The observed 
correlation pattern in central collisions is quite different 
from the expectations for statistical global momentum 
conservation \2]\ . Figure EJb) shows the diagonal and 
off-diagonal projections of the away-side three-particle 
correlation result from the 12% central data. The off- 
diagonal projection of our result with a = 6 = 1 is also 
shown. The off-diagonal side peaks are prominent; these 
peaks are evidence of conical emission of charged hadrons 
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correlated with high p± trigger particles. The side peaks 
in the diagonal projection contain conical emission and 
possibly other contributions such as k± broadening, large 
angle gluon radiation, and deflected jets. 

The angular distance 9 of the off-diagonal peak loca- 
tions from TT is obtained by fltting the off-diagonal projec- 
tions to a central plus two symmetric side Gaussians. For 
12% central Au-|-Au, 6l=1.37±0.02(stat)t5J;[!^(syst) ra- 
dian. The difference between 6 and the fit position 
to two-particle correlation may arise because the lat- 
ter measures a combination of effects. The value 
of 6 does not depend on centrality or the associ- 
ated particle p±. For p±— 0.5-1, 1-1.5, 1.5-2, 2- 
3 GeV/c, 6i=1.38±0.03(stat)l^{!^(syst), 1.36±0.04^°:^^, 



1.29±0.04ro:ro' and 1.31±0.05I^ o^, respectively. If the 
observed conical emission is generated by Mach-cone 
shock waves, the measured angle 9 reflects the speed of 
sound in the created medium averaged over the evolution 
of the collision M, M . 
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FIG. 4: Three-particle correlation strength per radian^ versus 
where A^part is the number of participants. Some 
of the data points have been displaced in A'part for clarity. 

To characterize the correlation strength, the aver- 
age signals are evaluated within 0.7x0.7 radian^ cen- 
tered at (A(/>i,A02) = (7r,7r), (vr ± 1.37, tt ± 1.37), and 
(tt ± 1.37, TT =F 1.37). Figure 2] shows the average signal 
strength 19] in pp, d-|-Au, and Au-|-Au collisions as a 
function of (A'part /2)^^'^. The signal strength increases 
and appears to saturate in central collisions. While the 
away central peak is the dominant structure in pp, d-|- Au, 
and peripheral Au-f Au, the diagonal and off-diagonal 
side peaks increase rapidly in strength with centrality 
and become the dominant structures in central Au-f Au 
collisions. 

In conclusion, the flrst three-particle azimuthal cor- 
relation measurements with a high transverse momen- 
tum trigger particle are reported by the STAR exper- 
iment. The analysis treats the event as composed of 
two components, one correlated with the trigger and the 
other, background. Results are presented for minimum 
bias pp, d-f Au, and different centralities in Au-|-Au col- 
lisions at y/Sj^^ — 200 GeV between a trigger particle 



of 3 < < 4 GeV/c and two associated particles of 
1 < < 2 GeV/c. Dijet structures are observed in 
pp, d-|-Au and peripheral Au+Au collisions, with a pro- 
gressive diagonal elogation of the away-side central peak. 
Distinct peaks at 0=1.37 ± 0.02(stat)^p p^(syst) from tt 
are observed on the away side in central Au-|-Au colli- 
sions, with correlated hadron pairs far apart, symmetric 
about TT. These structures are evidence of conical emis- 
sion of hadrons correlated with high p± particles. The 
conical emission angle is measured to be independent of 
the associated particle p± . 
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